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Abstract. The production of the stable isotope ^Li in standard Big Bang nucleosynthesis has recently 
attracted much interest. Recent observations in metal-poor stars suggest that a cosmological ®Li plateau 
may exist. If true, this plateau would come in addition to the well-known Spite plateau of ^Li abundances 
and would point to a predominantly primordial origin of ^Li, contrary to the results of standard Big Bang 
nucleosynthesis calculations. Therefore, the nuclear physics underlying Big Bang ^Li production must be 
revisited. The main production channel for ^Li in the Big Bang is the ^H(a,7)^Li reaction. The present 
work reports on neutron-induced effects in a high-purity germanium detector that were encountered in 
a new study of this reaction. In the experiment, an a-beam from the underground accelerator LUNA 
in Gran Sasso, Italy, and a windowless deuterium gas target are used. A low neutron flux is induced by 
energetic deuterons from elastic scattering and, subsequently, the ^H(d,n)^He reaction. Due to the ultra-low 
laboratory neutron background at LUNA, the effect of this weak flux of 2-3 MeV neutrons on well-shielded 
high-purity germanium detectors has been studied in detail. Data have been taken at 280 and 400 keV 
a-beam energy and for comparison also using an americium-bery Ilium neutron source. 
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1 Introduction 

1.1 Cosmology and ^Li 

For 30 years now, the cosmological lithium puzzle has frus- 
trated the efforts of observers and cosmologists [l|, a re- 
cent review] . Observations in metal-poor stars consistently 
show a factor of two or three lower abundance of the main 
stable lithium isotope ^Li (isotopic abundance on Earth 
92.41%) than what is predicted by standard Big Bang nu- 
cleosynthesis, even though some possible stellar solutions 
have been suggested [1]. 
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Recent observations of the second stable lithium iso- 
tope ^Li (isotopic abundance on Earth 7.59%) in metal- 
poor stars have introduced a possible additional puzzle, 
this time concerning cosmological ^Li. The observed 
^Li/^Li abundance ratio of about 0.05 [1] largely exceeds 
the standard Big Bang nucleosynthesis prediction. Even 
though many of the claimed ^Li detections may be in er- 
ror [3, Q , for a few metal-poor stars @ there still seems 
to be a ^Li isotopic abundance of a few percent Q • These 
observations are much higher than the predicted ^Li yield 
from standard Big Bang nucleosynthesis Q. 

^Li is much more easily depleted through nuclear re- 
actions than its more abundant sister isotope ^Li 0, and 
non-cosmological scenarios for significant ^Li production 
all co-produce ''Li and thus worsen the ''Li problem [T^, 
a very recent example]. Therefore, the detections of ^Li 
in very old stars raise the question of a possible cosmo- 
logical production. However, with standard Big Bang nu- 
cleosynthesis producing much too little ^Li Q, all cosmo- 
logical scenarios yielding enough ^Li involve non-standard 
physics plI - IT^ . Before studying such exotic scenarios, it 
is important to first determine precisely how much ^Li can 
be produced in standard Big Bang nucleosynthesis. With 
this experimental baseline, any missing additional ^Li pro- 
vided by non-standard approaches can be quantified. 

1.2 The nuclear physics of ^Li production 

The ^H(a,7)^Li reaction is the dominant nuclear reaction 
for ^Li production in standard Big Bang nucleosynthesis 
0. At the energicfQ relevant to Big Bang nucleosynthesis, 
E w 50-300 keV, the ^H(a,7)^Li cross section is very small. 
Therefore, it has never been measured experimentally at 
such low energies, and theoretical predictions remain un- 
certain |15l - ll7l |. In the low-energy domain, the cross sec- 
tion (T2a{E) is usually parameterized as the astrophysical 
S-factor S2a{E) defined by 

S2i{E) =Ea2i{E) exp(27r7,(£;)) (1) 

with E the center-of-mass energy, and 2TTrj{E) the Som- 
merfeld parameter parameterizing the exponential-like en- 
ergy dependence of the probability of tunneling through 
the Coulomb barrier [Tl] . 

The reaction has been studied previously by in-beam 
7-spectrometry around the E = 0.711 MeV resonance [l^- 
At even higher energies, there are data by in-beam de- 
tection of the ^Li reaction products [2^. In both cases 
[TqI . [20| , an a-beam had been incident on a deuterium gas 
target. A third experiment at very low energy, using a 
deuterated polyethylene target, just resulted in an upper 

limit mi. 

Two attempts to determine the ^II(a,7)^Li cross sec- 
tion have been made using the Coulomb dissociation tech- 
nique [l^, [12], in which an energetic ^Li beam is shot 
on a target of high nuclear charge. The time-reversed re- 
action ^Li(7,Q;)^H is then studied using virtual photons. 

^ Throughout this text, E always denotes the center-of-mass 
energy, and Ea the a-beam energy in the laboratory. 



The cross section of interest is obtained by applying the 
detailed balance theorem [2^. This method is especially 
sensitive to quadrupole (E2) transitions in the excited ^Li 
nucleus, largely neglecting dipole transitions. It is further- 
more limited by possible background from non-Coulomb, 
i.e. nuclear breakup 

Such experiments have been carried out recently, us- 
ing a 150MeV/A ^Li beam and the KaoS spectrometer at 
GSI [3 • It was found that nuclear breakup dominated the 
signal, so only an experimental upper limit could be de- 
rived |16| . but the observed angular distributions seemed 
to support a theoretical excitation function developed in 
the same paper. It is stated that the previous Coulomb 
breakup data [22[, which had been obtained at a much 
lower ^Li projectile energy of 25MeV/A, are probably 
even more strongly affected by nuclear breakup, so that 
also these data should be interpreted as upper limits [l6| . 

Therefore, direct data on the ^H(a,7)^Li cross section 
at Big Bang energies are still needed. The present work 
lays the foundation for such an experiment at the Lab- 
oratory for Underground Nuclear Astrophysics (LUNA, 
mill) 400 kV accelerator. 



1.3 Processes expected to take place in the target 

In order to prepare the planned ^H(a,7)^Li experiment, 
it is necessary to study in details the dominating back- 
ground. Inside the target, mainly the following nuclear 
reactions are expected to take place: 

2H(a, -ffU Q = 1.474 MeV (Rl) 

2H(a,a)2H — > ^}l{d,nfm g = 3.269 MeV (R2) 
— > p)3jj g = 4.033 MeV (R3) 

Along with the reactions, also their g-values are listed. 
Reaction Rl is the main reaction to be studied. Reac- 
tions R2 and R3 are parasitic reactions induced by elas- 
tically scattered deuterons. As these deuterons have a 
maximum energy of only 356 keV (for E^ = 400 keV), 
other deuteron-induced reactions are negligible due to sup- 
pression by the Coulomb barrier. However, there is no 
Coulomb barrier for the energetic neutrons released by 
reaction R2, and the protons released by reaction R3 have 
an average energy of 3 MeV. 

Therefore, these protons and neutrons may give rise 
to further reactions on the structural material of the gas 
target system and on the detector material. The cross sec- 
tions of these two conjugate reactions are similar, and they 
have nearly the same energy dependence at the energies 
relevant here [13]. As a consequence, the easily detectable 
protons from reaction R3 may in principle be used to mon- 
itor, in a relative manner, also the yield of reaction R2. 



1.4 Neutron-induced effects in germanium detectors 

There is a considerable body of literature on this prob- 
lem. Soon after the advent of lithium-drifted germanium. 
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or Gc(Li), detectors, the effects of neutrons on these de- 
vices has been studied [H, HI]. The energy loss of re- 
coiling germanium atoms in the germanium explains the 
characteristic tr iang ular shape due to (n,n'7) reactions in 
germanium [s^, l3lj | These effects have already been in- 
cluded in Monte Carlo simulations [111, and the possi- 
ble discrimination against unwanted neutron effects in a 
highly segmented detector array has been studied [33| . 
The response of HPGe detectors to strong Am-Be neu- 
tron sources [S^, HI] and the effects of cosmic-ray induced 
neutrons on low-background germanium detector systems 
[36l - [38l ] have both been studied in detail before. 

In the present work, a detailed study of the effects 
of low neutron fluxes on an underground, well-shielded 
germanium detector is presented, using parasitic neutrons 
produced by the ^H(d,n)'^He reaction in the in-beam ex- 
periment and, subsequently, an americium-beryllium (Am- 
Be) neutron source. 

This work is organized as follows. The experimental 
setup is described in secjU the Monte Carlo simulations 
and their validation using an americium-beryllium neu- 
tron source in sec.[3l The main observations with the sili- 
con and the germanium detectors, detecting charged par- 
ticles and 7 rays, respectively, during the in-beam exper- 
iments are summarized in sees. HI and \5\ The results and 
the effects for the ^H(Q;,7)^Li experiment at LUNA are 
discussed in sec.|6l and a summary is given in sec.[7l 



2 Experimental setup 

2.1 The LUNA 400 kV underground accelerator 

The experiment was performed at the Laboratory for Un- 
derground Nuclear Astrophysics (LUNA) 400 kV accelera- 
tor in the Gran Sasso laboratory, Italy (LNGS). The Gran 
Sasso underground facility is shielded from cosmic-ray in- 
duced muons by the Gran Sasso mountain, reducing the 
muon flux by six orders of magnitude [39l ]. The ambient 
flux of energetic neutrons at Gran Sasso is of the order 
of 10~^s~^ cm~^ 113, nil, three orders of magnitude lower 
than at the surface of the Earth. 

LUNA is the world's only underground accelerator fa- 
cility. It is dedicated to cross section measurements of as- 
trophysicall y re levant reactions well below the Coulomb 
barrier 25. [2611 . taking advantage of its ultra- low back- 
ground 42h44|. 

During the experiment described here, the LUNA 
400 kV accelerator [i^ provided a *He^ beam of 280 and 



400 keV energy, with an intensity of 100-350 /iA (typical 
value 300/^ A). 



2.2 Windowless gas target system 

After the analyzing magnet, the ^He"'" beam was directed 
to a windowless gas targ et sy stem consisting of three dif- 
ferential pumping stages |46l ] . After passing through a se- 
ries of three collimators of 25, 15, and 7 mm diameter, the 



beam enters the target chamber (fig.P). The final collima- 
tor with 7 mm diameter is directly at the entrance of the 
target chamber. The other two collimators are upstream 
at 69 and 51cm distance, respectively. Inside the target, 
deuterium gas of 99.8% chemical and isotopic purity was 
maintained at a constant pressure of O.Bmbar via a capac- 
itive pressure sensor (independent of gas type, calibrated 
to 0.25% precision) controlling the gas inlet via a feedback 
loop. 

The target chamber is a box-shaped container of 44 cm 
length and 11x13 cm area, made of 2 mm thick AISI 304 
steel. 6 cm downstream of the initial wall, a recess of 13 cm 
length and 4.4 cm height is introduced in order to insert 
the HPGe detector close to the beam line. At the center of 
the target chamber, a smaller, box-shaped AISI 304 steel 
inner tube of 1.8 x 1.8 cm area and 17 cm length with a wall 
thickness of 1 mm was introduced. This inner tube lim- 
its the lateral length that elastically scattered deuterons 
travel inside the gas, and therefore the deuterium target 
for the ^H(d,n)^He reaction. The inner tube and the inner 
part of the collimator at the entrance of the target cham- 
ber were made from interchangeable parts, so they were 
replaced with a fresh component from time to time, typ- 
ically after about ~50G C of beam charge. Due to proper 
inner linings, the position of the components was repro- 
ducible to better than 1 mm. 

Outside the target, the differential pumping system 
maintains already a medium vacuum (lO^'^mbar) at the 
first pumping stage, evacuated by a 2050 /h Roots pump. 
The second pumping stage is equipped with three 1500 1/s 
turbomolecular pumps, and the third pumping stage with 
a 360 1/s turbomolecular pump. The vacuum in the sec- 
ond and third pumping stages is in the 10~^-10~^ mbar 
range. The backing pumps for the three pumping stages 
are oil-free rotary vane compressors. 

After traveling for 17.7 cm through the gas target, the 
ion beam is stopped on the copper head of a beam calorime- 
ter with constant temperature gradient to measure the 
beam intensity . Resistive heaters controlled via a feed- 
back loop and forced cooling maintain a constant temper- 
ature gradient between 70 °C on the hot side, facing the 
beam, and 0°C on the cold side of the calorimeter. 

The gas temperature inside the target chamber varied 
monotonously between 16 °C next to the water-cooled en- 
trance collimator and 70 °C immediately adjacent to the 
hot side of the calorimeter. The pressure profile had been 
measured in a similar setup in a previous experiment at 
LUNA and was found to be flat [43]. It is assumed here 
that it is again flat. 

The effective target density depends on target pressure 
and temperature, but also on the beam heating effect [48]. 
This effect has been studied previously at LUNA for ^He 
gas, in a setup similar to the present one, by double elastic 
scattering [4§|- For ^He gas, the elastic scattering data 
[49| are in agreement with the predicted beam heating 
correction AT from the beam current /, energy loss in the 
target dE/dx, heat conductivity A and the approximate 
radii of target chamber rchambcr and beam rBcam using 
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Fig. 1. Experimental setup, as seen from the side. The central chamber of the windowless gas target is seen near the center 
of the plot. The a-beam enters the target from the left through a 4 cm long collimator of 7 mm inner diameter and is stopped 
on a massive copper beam calorimeter. The germanium (for 7 rays) and silicon (for charged particles) detectors are also shown, 
as well as the PtlOO temperature sensors and the tube leading to the capacitive pressure sensor. The setup is surrounded by a 
lead castle and walls of borated polyethylene. The inner lead castle is surrounded by an anti-radon box made of acrylic glass. 



the following formula [5C 



AT 



dE 



dx 2ttX 



1 1 ^Chamber 

In 



(2) 



Using / = 200-300 yuA, dE/dx = 0.17 (0.19) keV/cm at 400 
(280) keV, A = 1.31 niW/(Kcni), rchambcr « 9nim, rBcam 
« 3.5 mm, a correction of 4-7 K is obtained, correspond- 
ing to a 1-2% reduction in effective target density. As in 
previous work [49l |. a conservative relative uncertainty of 
40% is assumed for this correction. 



2.3 Silicon particle detector 

For the detection of protons from reaction R3, a 1500 /im 
thick silicon detector of 450 mm^ active area is installed 
atop the center of the gas target chamber (fig. HI . At this 
position, it can detect protons from reaction R3 taking 
place both within the gas volume and in the wall of the in- 
ner tube, within the limits of the detector solid angle. The 



silicon detector is covered with a 25 /^to thick aluminum 
foil in order to stop elastically scattered ^He and parti- 
cles (with a range 1.5 and 3 /xm, respectively) while letting 
the energetic protons pass. Without the foil, the measured 
energy resolution of the detector itself is about 25 keV full 
width at half maximum (FWHM) for 5.5 MeV a particles 
from an ^'^"'^Am source. 

Without active cooling and at high beam intensities, 
the silicon detector reaches a temperature of up to 45 °C, 
as measured by a PtlOO resistor at the detector enclo- 
sure. This is so, because convective cooling is negligible at 
the present low target gas pressure. This increase in tem- 
perature leads to degraded performance of the detector, 
because the main contributor to its energy resolution is 
thermal noise. 

Therefore, an active cooling is provided by a Peltier 
element on the back side of the detector. Its hot side is 
thermally connected by an aluminum bar to the target 
chamber wall (fig. [I]). When cooled by the Peltier element, 
the silicon detector temperature never exceeds 23 °C, the 
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recommended working temperature specified in the detec- 
tor data sheet. 



2.4 Germanium 7-ray detector and lead shield 

The 7 rays are detected in a large high-purity germanium 
detector of 137% relative efficiency. This detector is placed 
inside the recess in the target chamber (fig.[l]). The sym- 
metry axis of the HPGe detector is 7.5 cm downstream of 
the entrance of the gas target, and its end cap just 1.5 cm 
below the beam axis. This detector has been optimized 
for ultra-low background, and it has previously been used 
to construct a setup with the lowest background ever ob- 
served for an in-beam 7-spectrometry setup [43l |. 

The 7-ray detection efhciency was measured with cal- 
ibrated -'^•^^Cs, ^°Co, and ^'^Y sources (figH]) as a function 
of their position along the beam axis. Due to the close 
geometry used, the data had to be corrected for the true- 
coincidence summing out effect, a correction amounting 
to up to 12% at the closest distance studied here. 

The typical resolution observed for the 1.333 MeV line 
of ^"^Co is 2.7 keV when the gas target pumps were switched 
off, a few tenths of a keV worse when the pumps were run- 
ning. For the same 7-line, with the source placed near the 
center of the detector [x = 7.5 cm) a peak to Compton 
ratio of 65 was found. 

The detector was shielded from environmental radionu- 
clides by a lead shield of at least 20 cm thickness (fig.[I|). 
The lead is selected for a ^^°Pb content of just 25 Bq/kg 
(43| . In order to limit the impact of the neutrons, the setup 
was enclosed by a 5 cm thick shield of borated polyethy- 
lene (fig.[l]), in order to exclude any possible background 
for other experiments in the Gran Sasso underground fa- 
cility. This shield attenuated the additional neutron back- 
ground produced by the present experiment outside the 
LUNA experimental hall to a level that was well below 
the already ultra-low ambient neutron flux at Gran Sasso 
[40l . l4l| . In order to reduce the background from radon 
decay products, the setup was enclosed in an anti-radon 
box made of acrylic glass that is flushed with radon-free 
nitrogen evaporated from a Dewar vessel. 



2.5 Data acquisition 

The data from the germanium detector were passed on to 
two independent data acquisition systems for processing 
and storage. One branch of the data included the TNT2 
digital data acquisition system in a CAEN N1728B mod- 
ule, operated in list mode. The dead time was estimated 
using a pulser. 

In the second branch, histograms were recorded at 
regular intervals using the Ortec 919E EtherNIM ana- 
log to digital converter and multichannel buffer. Here, the 
Gedckc-Hale algorithm was used for dead time estimation. 
The data from the silicon detector were treated just by the 
Ortec branch, with a dead time of less than 1% estimated 
with a pulser connected to the preamplifier test input. 




Data, 662 keV (Cs-137) 
Data, 1333 keV (Co-60) 
Data, 1836 keV (Y-88) 
Monte Carlo, 662 keV 
Monte Carlo, 1333 keV 
Monte Carlo, 1836 keV 



5 10 15 

Position after end of collimator, x [cm] 

Fig. 2. Full-energy peak detection efficiency measured with 
several calibrated 7-ray sources, as a function of the position 
X after the end of the final collimator. The efficiency from the 
GEANT4 Monte Carlo simulation is also shown for several 
positions. 



The main parameters of the gas target system and 
calorimeter, such as the pressures observed in the various 
pumping stages and the temperatures at several places 
inside the calorimeter, were logged at intervals of a few 
seconds via a Labview-based slow control system. 



3 Monte Carlo simulation 

The experimental setup was coded in two versions of 
GEANT, each of them suited to the special needs of the 
case: A specially adapted version of GEANT3 taking into 
account the correct energy loss for charged particles was 
used to model the data in the silicon particle detectors. 
The details of this code have been discussed previously 
[sH ] , and the output from the GEANT3 simulation is com- 
pared with the data in sec. SI 

The response to 7 rays and neutrons was coded in 
GEANT4 j5j|. The GEANT4 simulation is described in 
the present section. 



3.1 GEANT4 simulation 

The GEANT4 version used was 4.9.2, and the 
QGSP_BIC_HP physics list was activated. This physics list 
uses the GEANT4 binary cascade (BIG) mechanism for 
primary protons and neutrons with energies below 10 GeV. 
The excited nucleus created by the primary nucleon is then 
passed to the precompound model, which describes the nu- 
clear de-excitation. For neutrons below 20 MeV down to 
thermal energies, the data-driven high-precision neutron 
package (NeutronHP) is used. The NeutronHP version used 
was G4NDL3.13, which is based on the ENDF/B-VI and 
JENDL libraries. The high-energy interactions are taken 
from the JENDL High Energy File 2007. 

The first step of the GEANT4 simulation was to ad- 
just the details of the geometry for the HPGe detector 
used here, starting from the nominal geometry, until the 



6 M. Anders et al. (LUNA coUab.): Neutron-induced background by an a-beam incident on a deuterium gas target. 



data from the radioactive 7-ray sources ^^^Cs, ^°Co, and 
placed at various positions along the beam path were 
reproduced to a precision of better than 10%. The crystal 
diameter and length were reduced after a repair and etch- 
ing process of the detector used here. In order to reproduce 
the measured efficiency profile (figH]), the following values 
were applied, with the data sheet values for the detector 
before repair given in brackets: Crystal diameter 86 mm 
(91mm), length 85mm (91mm), distance to the casing 
12 mm (4 mm). The quoted relative efficiency of the de- 
tector before (after) repair was 150% (137%). With the 
adopted dimensions given here, the simulation predicts a 
relative efficiency of 130%, close to the data sheet effi- 
ciency after repair. 

For a ^"Co source, a peak to Compton ratio of 64 was 
found in the simulation near the center of the target (x 
= 7.5 cm, with x the position measured after the end of 
the final collimator), after scaling to obtain the same res- 
olution as in the experiment. This is consistent with the 
experimental value of 65 (sec. l2.4p . showing that the sim- 
ulated geometry seems to be appropriate. Also at other 
positions and for other sources, the efficiency data with 
7-ray sources was well reproduced by the Monte Carlo 
simulation (fig HI • 

As a second step, the neutron-induced background in 
the in-beam experiment was modeled in the given geom- 
etry. This was done in several discrete steps, in order to 
limit the computation time: 

1. *He+ ions of the correct laboratory energy (280 or 400 
keV) are shot into the gas target, and the Rutherford 
scattering reaction ^H(a,a)^H is implemented. The en- 
ergy spectrum of emitted deuterons, integrated over all 
angles, is recorded and kept as input for the following 
step. 

2. A source of energetic deuterons, distributed in energy 
according to their previously derived energy spectrum, 
is created. The point of emission is a cylinder of the 
dimensions of the He^ beam, along the beam axis be- 
tween final collimator and beam stop. The radial pro- 
file of emission points of origin is randomly distributed 
with a Gaussian distribution following the beam pro- 
file, and the angle of emission is also randomly dis- 
tributed. The deuterium gas target is then irradiated 
with this source, and the neutron-emitting reaction 
^H(d,n)'^He is implemented. The interaction of the re- 
sulting neutrons with the structural and shielding ma- 
terials and with the germanium detector material is 
then simulated. 

3. The energy deposited in the HPGe detector by inci- 
dent 7 rays is then recorded and histogrammed. The 
energy deposited by recoiling germanium atoms inside 
the detector crystal is also included, taking into ac- 
count quenching effects [3^ . 

3.2 Benchmarking of the GEANT4 simulation using an 
americium-beryllium neutron source 

In order to benchmark the GEANT4 simulation, a weak 
americium-beryllium (Am-Be) neutron source, emitting 
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Fig. 3. Spectrum of the AmBe neutron source from the ex- 
periment and from the Monte Carlo simulation. 

MeV neutrons by the ^Be(Q!,n)^^C reaction, was intro- 
duced to the setup. The Am activity was 185 kBq, lead- 
ing to an estimated neutron source strength of 13n/s. The 
^^^Am nitrate was contained in a thin layer between two 
0.1mm thick beryllium disks of 11mm diameter, housed 
in a 5 mm thick polyethylene cylinder of 18 mm diameter. 

The Am-Be source was placed in the center of the tar- 
get chamber. Due to restrictions on the use of neutron 
sources in the LNGS underground facility, the running 
time was only 10 hours, limiting the statistical precision 
of the data. Therefore, the data had to be rebinned in or- 
der to make the spectra comparable, and any comparison 
by necessity concentrates on the general features of the 
spectrum (fig.[3|). 

The triangular feature at 691 keV is not well repro- 
duced in the simulation, for reasons discussed below (sec. (5]) 
At lower energies, < 500 keV, the simulation signifi- 
cantly underpredicts the data. This is due to the Compton 
continuum from the 7 rays at 662 and 722 keV that are 
emitted by the ^^^Am in the Am-Be source. These weak 
7-ray branches of ^^^Am arc not included in the Monte 
Carlo simulation. In addition to energetic neutrons from 
the ^Bc(Q!,n)^^C reaction, an Am-Be source also emits 
Doppler-broadened 4.4 MeV 7 rays from the decay of the 
first excited state of ^^C. These 7 rays were not included 
in the Monte Carlo simulation, but they have only a very 
limited influence near the ROIs for the ^II(a,7)^Li reac- 
tion. 

Overall, the spectrum from the Am-Be source is well 
reproduced by the simulation, including the neutron tri- 
angles except for the one at 691 keV where there seem to 
be problems with the treatment of the internal conversion 
process in GEANT4 (secO. 

4 In-beam observations with the silicon 
particle detector 

Typical spectra taken with the silicon detector are shown 
in flg.m for the two a-bcam energies used here. The spec- 
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Fig. 4. Experimental (red solid line) and GEANTS-simulated 
(blue dashed line) particle spectrum in the silicon detector. 
Top (bottom) panel, runs at 400 (280) keV. The simulated 
spectrum is given in arbitrary units. Due to the strong pream- 
plification, the spectra are affected by electronic noise below 
IMeV. 



tra show a broad peak at about 2.0 MeV, corresponding 
to the most probable energy of the protons emitted by 
reaction R3, after ?«0.5MeV energy loss in the 25 //m alu- 
minum foil. The width of the peak is given by the kinemat- 
ics and energy straggling in the Al foil. For Ea = 400 kcV, 
the high-energy cutoff is found at 3.0 MeV, as expected 
from geometry and kinematics. 

It is evident that the energy of the proton peak due 
to reaction R3 is almost independent of beam energy, as 
it is dominated by the reaction Q-value. Similar consider- 
ations apply to the neutrons from reaction R2, which are 
not directly detected here. The Monte Carlo simulation 
reproduces the shape of the experimental data fairly well 
(fig.®. 

It is expected that some of the energetic deuterons 
created in the beam are implanted in the walls of the inner 
tube, at a depth of up to their 1.6 /im range. There, they 
may again form a target for reactions R2 and R3. The 
protons created in nuclear reactions in this implantation 
layer will have slightly lower energy than protons created 
in the gas, due to the slowing of the deuterons in this layer 
before the reaction takes place, and due to the proton 
energy loss (up to 0.07 MeV) in the implantation layer. 
However, this effect is much smaller than the width of the 
peak. 

During the experiment, it was found that the proton 
rate measured with the silicon detector increased with ac- 
cumulated charge, by about 10% for each step of 100 C 
of charge. This confirms that the deuterium implanted in 



the target chamber actually does contribute to the proton 
yield, and by extension also the neutron yield. 



5 Observations with the germanium 7-ray 
detector with and without beam 

5.1 Adopted 7-ray regions of interest (ROIs) 

For the planned study of the ^II(Q;,7)^Li reaction with its 
very low expected yield, it is foreseen to use the highest 
beam energy available at LUNA, Ea ~ 400 keV. For sys- 
tematic checks, a second, lower beam energy is selected at 
280 keV, chosen so that the regions of interest (ROIs) do 
not overlap for the two reactions. 

The ROIs are determined as follows. The ^II(a,7)^Li 
reaction gives rise to a single 7 ray, that is observed in the 
detector at energy E^ 



E^ = Q + E±AE-_ 



Dopplcr 



AEi 



Recoil 



(3) 



The 7-ray energy shift due to the recoiling compound nu- 
cleus is negligible here, Z\i?ptocoii = 0.2 keV. The Doppler 
correction, however, is significant, with the full Doppler 
shift amounting to ziE'Doppier ~ 16keV at E^ = 400 keV. 
The energy loss over the length of the target amounts to 
about 2keV. 

As the target is extended over and beyond the full di- 
ameter of the germanium detector leading to emissions 
before and behind the detector (fig. [J), the 7 rays, includ- 
ing their flanks, fall into a ROI that ranges from 1552.5- 
1581. 5keV for E^ = 280 keV, and from 1589.5-1624.3keV 
for Ea ~ 400 keV. These ROIs will be assumed throughout 
the present work. 

The respective contributions of electric dipole and elec- 
tric quadrupole capture to the cross section are known 
only from theory (Tsl . IT6| . Therefore, the angular distribu- 
tion of the emitted 7 rays is highly uncertain. Because of 
the strong Doppler effect discussed above, the unknown 
angular distribution translates into an unknown shape of 
the 7 peak. In order to be insensitive to the transition 
type (dipole or quadrupole), the entire region of interest 
is used here. 



5.2 Laboratory background 

Laboratory background induced by cosmic rays is so low 
inside the deep- underground LUNA facility that it is negli- 
gible for the purposes of the present study. However, this is 
not the case for background from environmental radionu- 
chdes, such as '^^K and the nuchdes of the ^^®U and ^^^Th 
decay chain. This latter background was attenuated by 
the lead shield (sec. l2.4p . and the laboratory background 
counting rate for 200 keV < E^ < 2700 keV was found to 
be reduced by a factor of 800 compared to the unshielded 
setup. Inside the ^H(Q;,7)^Li region of interest, the reduc- 
tion reaches a factor of 1100. 

It has to be noted that a previous version of the present 
setup had achieved an even lower background level, mainly 
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Fig. 5. Spectra taken with the germanium detector. Blue full line: in-beam spectrum at Ea ~ 400 keV, pTargct = O.Smbar, 
laboratory background subtracted. The quantity plotted is the counting rate [h^^ keV^^]. In order to obtain the yield [C~^ 
keV"'^], the values plotted have to be divided by 1.12 C/h. Grey thin line; Laboratory background. The most important in-beam 
lines due to (n,n'7) and (n,7) processes on structural and shielding materials are marked with arrows, and the relevant target 
nuclide is given, as well as the 7-ray energy in keV. This information is also summarized in table [T] The asymmetric lines due 
to (n,n'7) processes in the germanium detector itself are listed in table O 



by lining the inside of the shield with 4 cm of eleetrolytie 
eopper, attenuating 7 rays and bremsstrahlung from the 
. For the present purposes, however. 



remaining ^Pb 
no copper liner was used, because of the high neutron 
capture cross section of copper. The lesser shielding with 
respect to Ref. [1^ leads to a factor of 16 higher laboratory 
background in the ROI, but it will be shown below that 
this no-beam background is still acceptable for the present 
purposes. 

Two 7 rays in the remaining laboratory background 
fall inside the region of interest: The Ac line at 1588 keV 
and the ^1231 line at 1621 keV, both are part of the '^^'^Th 
decay chain. The ^^^^Th chain nuclides are believed to be 
an external contamination of the lead used for the shield- 
ing e.g. by dust particles. 

At higher 7-ray energies, the 2614 keV line of ^os^pj jg 
clearly visible: again a member of the ^'^^Th decay chain. 
This line can in principle display a double escape peak 
inside the region of interest. However, the large size of the 
detector ensures that double escape peaks are negligible 
for all practical purposes. Even still, the Compton con- 
tinuum caused by this line clearly influences the counting 
rate in the region of interest. The same is true for several 
lines from the ^^^Rn daughter ^^'^Bi. in particular at 1730 
and 1764keV. 

Due to a constant flushing of the anti-radon box with 
dry nitrogen, the laboratory background in the two 



'^ll{a,j)^Li ROIs was found to be stable within statistics 
over a period of several months, at a level of 
(2.40±0.16) X 10-2 and (3.32±0.17) x lO'^ counts keV'^ h"' 
for the Ea = 280 and 400 kcV ROL respectively. 

In conclusion, the level of laboratory background is 
sufhcicntly low for the present purposes. It is stable in time 
and does not depend on the beam, so it can be subtracted. 



5.3 Structure of the background induced by the ion 
beam 

Due to the parasitic reactions R2 and R3, a certain amount 
of unwanted energetic neutrons, protons, tritons, and "^He 
particles are produced. They, in turn, may be captured by 
some structural, shielding or detector material. 

Energetic neutrons cause two main categories of effects 
in the setup. The first consists of capture and scattering ef- 
fects on structural and shielding materials. They give rise 
to a number of well-defined 7 lines with Gaussian shape 
(fig. O and table[T|). The detector end-cap consists of elec- 
trolytic copper, and the beam stop also consists of copper. 
A number of lines from the two stable copper isotopes 
63,65q^ are indeed observed (red markers in fig. [5]) |53l |. 
The gas target chamber, the collimator and the small tube 
inserted in the chamber consist of AISI 304 steel, and re- 
lated neutron-induced lines from scattering on 54,56pg j5g[ | 
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Table 1. 7-ray energies, in keV, of observed Gaussian-shaped 
neutron- induced lines (arrows in fig.[5ll , and the transition the 
lines originate from. See table[2] for the broadened Ge(n,n'7) 
peaks, and see the text for a discussion of the 198 keV line. 



^®Ni and ^^Cr are apparent (grey and turquoise markers 
in fig. [5]). At 803 keV, a line from neutron scattering on 
206p|^ from the massive lead shield surrounding the target 
ehamber is observed. 

The second main effect is caused by (n,n'7) processes 
of energetic neutrons on the germanium detector material 
itself. In addition to the emission of a 7 ray from the dc- 
excitation of an excited nuclear state in the germanium 
target nucleus, the recoiling tar get nucleus also deposits 
energy in the detector material [32|, leading to a charac- 
teristic triangular shape starting at the energy of the 7 
ray. The broad features due to this effect (fig.E]) can be 
attributed to a number of excited states in several sta- 
ble germanium isotopes (table [2]) . A number of the fea- 
tures observed here have been reported in previous work 

The ''^™Ge metastable state at = 198 keV is a spe - 
cial case. It is excited by the ™Ge(n,7) reaction pol. [BGf. 
Due to its long half-life of 20 ms, ^^™Ge decays only after 
the recoiling nucleus has been stopped. The ^^™Ge level 
at 198 keV decays via the 175 keV state in ^^Ge and on- 
ward to the ground state. The sum of these two 7 rays 
gives rise to a Gaussian peak at 198 keV 7-ray energy in 
the spectrum (fig. [5]). 

In principle, further nuclear reactions in addition to 
the neutron-induced reactions discussed here are possible. 
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Table 2. 7-ray energies, in keV, of broadened Ge(n,n'7) peaks 
in the 200-2200 keV energy range (fig.[S|. The germanium nu- 
clear level energies are taken from the ENSDF file [H^. Pre- 
vious observations of these lines are indicated, where they are 
known to the authors. The 691 keV feature is caused by internal 
conversion (IC). See table[T]for the Gaussian-shaped neutron- 
induced lines. 



In addition to energetic neutrons, reactions R2 and R3 
give rise to 0.8-1.0 MeV tritons and "^He particles. These 
particles, in turn, may in principle give rise to the 
^II(t,n)'^He and the ^H('^IIe,p)^He reactions when pass- 
ing the deuterium gas of the target, and to further triton- 
and "^Hc- induced reactions in the structural, shielding, and 
detector materials. Also (p,p'7) reactions are in princi- 
ple possible. However, the cross sections for these charged 
particle induced reactions are lower than those for neu- 
tron scattering, and no clear signature has been found for 
them in the 7-spectra, so they are neglected here. 

5.4 Overall energy dependence of the ion-beam 
Induced background 

The in-beam spectra display very similar features for the 
two beam energies Ea = 280 and 400 keV. The reason for 
this is that the maximum neutron energy varies not very 
much with beam energy, from 3.3 MeV at Ea = 400 keV 
to 3.1 MeV at Ea = 280 keV. The mutual similarity of 
the two neutron-induced spectra is actually even closer 
than the similarity between the experimental spectrum at 
Ea = 280 keV and the simulated spectrum at the same en- 
ergy (fig. [7]). Based on this information, it seems useful to 
attempt an empirical parameterization of the energy de- 
pendence of the ratio between the counting rates at these 
two beam energies, after normalization for beam intensity. 
If successful, such a parameterization will in principle al- 
low to use one of the two runs as a monitor run for the 
background of the other run, and vice versa, if the ROIs 
are disjunct. 

The ion beam induced 7 rays discussed in the previous 
section all give rise to a Compton continuum. Near the 
^II(a,7)^Li region of interest, this continuum is an order 
of magnitude higher than the laboratory background. 
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Fig. 6. In-beam spectra at Ea = 400 keV (blue line, laboratory background subtracted), Ea = 280 keV (green line, laboratory 
background subtracted), and laboratory background (black line). The flat regions used for the empirical parameterization 
fsec. l5.4l fig. [HI are marked by light blue boxes. The ROIs for Ea = 280 keV (400 keV) are indicated by red (blue) lines. For 
the in-beam data, the quantity plotted is the yield [C~^ keV~^]. In order to obtain the counting rate [h~^ keV~^], the values 
plotted have to be multiplied by 0.97 C/h (0.89 C/h) for Ea = 280 (400) keV. For the laboratory background, the counting rate 
[h"^ keV"^] is plotted. 



In the experimental spectrum at the two beam ener- 
gies considered here, flat regions have been selected which 
have a width of at least 15keV and do not exhibit any 
significant structure (fig-E]). In these flat regions, the spec- 
tral shape is nearly identical for the two cases Ea — 280 
and 400 keV, and the ratio of the two respective yields is 
near unity. However, for higher beam energy also a slightly 
harder neutron spectrum is expected, so the yield ratio 
should exhibit an energy dependence. This is indeed the 
case (fig-IH])- As expected, with increasing 7-ray energy the 
background yield at Ea = 400 keV increases slightly with 
respect to the 280 kcV yield. 



This behavior has been parameterized empirically with 
a quadratic function (fig.|Hl)- The quadratic function has 
been selected, because its x^/'^ value is 1.7, better than 
for a linear function which was also tried and which yields 
X^/i/ = 4.0. More complex fit functions have not been 



attempted here. 



6 Discussion 

In the following text, the in-beam data and the simula- 
tion are compared, and the consequences for the planned 
^II(a,7)^Li experiment at LUNA are discussed. 



6.1 Comparison of in-beam data and simulation 

The in-beam spectrum is plotted together with a Monte 
Carlo simulation for the Ea = 280 keV case (fig. [7]). As the 
spectra at Ea = 280 and 400 keV are very similar to each 
other (fig. HI), the following discussion also applies to the 
latter beam energy. 

The general features of the spectrum are rather well 
reproduced. Therefore, here only some notable but rel- 
atively minor discrepancies are discussed, following the 
numbering scheme in fig. [71 

1. The Gaussian sum peak denoting the ^^™Ge metastable 
state at E^ = 198 keV [29, Q&i is observed in the ex- 
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Fig. 7. Comparison of the experimental spectrum at Ea = 280 keV (green line) and the GEANT4 simulation for the same 
conditions (black line). The grey arrows are at the same positions as in fig.El but the labels have been omitted here for clarity. 
The numbered red arrows denote the most important discrepancies between the measured and the simulated spectrum and are 
discussed in the text. 



periment, but not in the simulation. This is due to 
the well known GEANT4 problem that while the de- 
cay of metastable states is correctly handled by the 
G4RDM module, they are not correctly produced by the 
inelastic capture module. The 198 keV state lives so 
long that the recoiling ''^™Ge is completely stopped 
long before it decays, so the typical neutron triangle 
structure does not apply. A peak at the same energy 
may also be caused by the ^^F(n,n'7) reaction, but this 
would require an unplausibly high amount of fluorine 
near the detector. 

The neutron triangle at 691 keV is much smaller in the 
simulation than in the data. The 0"*" excited state at 
i?x = 691 keV in ''^Ge decays exclusively by internal 
conversion to the 0+ ground state, meaning its energy 
is detected in the present large germanium detector 
with an efficiency close to 1. The Monte Carlo simula- 
tion erroneously treats this level as decaying by 7-ray 
emission, meaning the 7 ray is detected only with the 
intrinsic efficiency of the detector, which is consider- 
ably smaller than 1. 

The 1238 keV line assigned to the transition between 
the exited states at 2085 and 847 keV in ^^Fe shows 
up at 1233 keV instead. It should be noted that the 
decay of the 847 keV level to the ground state in ^^Fe 
is reproduced at the correct energy. The observed 5 keV 
7-ray energy difference is too large to be explained by 
the Doppler effect, which would amount to less than 
2keV in this case. This points to a possible error in 
the nuclear structure database included in GEANT4. 



4. There is a feature at 1632 keV in the simulation but 
not in the data. It looks like a Ge(n,n'7) triangle but 
with smaller width, thus lower deposited germanium 
energy. It is therefore most probably due to the pop- 
ulation of a relatively high-energy germanium level in 
the simulation but not in the data. 

5. Same as 4. above, but at 2240 keV. 

Owing to the general good agreement between data 
and simulation, no significant discrepancies can be found 
in the Compton continuum, especially not near the 
^H(a,7)^Li ROIs. This confirms that the simplified model 
introduced in sec. [3] correctly describes the main processes 
in the experiment. 



6.2 Effects for the data analysis of a ^H(a,7)^Li 
experiment 

If one assumes that the ^H(a,7)^Li cross section is equal 
to the recommended value from Rcf . [l^ , at Ea = 400 keV 
a signal to noise ratio of about 1:12 is obtained with the 
present neutron induced background. Therefore, it is im- 
portant that the background level is understood with a 
precision on the percent level. The present analysis may 
lay the groundwork for such a precision. 

It has been shown above that the Monte Carlo sim- 
ulation generally reproduces the experimentally observed 
spectrum rather well, both the general behavior and spe- 
cific features (fig.El). The ratio of the yields for the flat 
spectrum regions is in fair agreement with the data (fig. El) • 
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Fig. 8. Ratio of yields at Ea = 280 and 400 keV, respectively. 
The data are from the fiat regions marked in fig.IH] The em- 
pirical parameterization by a quadratic function is shown as 
a full line. The Monte-Carlo based ratios are scaled by a fac- 
tor of 1.02 in order to match the high statistics data point at 
200 keV. 

The mutual agreement in shape is even better between the 
Ea = 280 keV and 400 keV runs, respectively. 

One remaining discrepancy between the Ea = 280 and 
400 keV runs lies in the counting rate for the 1811 keV ^^Fe 
line. It is higher at 400 keV, possibly due to the some- 
what higher neutron energy for this run. The Compton 
edge for this particular 7 ray lies exactly between the two 
ROIs for the ^H(a,7)^Li reaction. Therefore, the effects 
of the Compton continuum and multi-Compton events of 
the 1811 keV 7 ray on these two ROIs has to be taken into 
account separately. Based on the data of the 1836 kcV 7 
ray from the source, it is estimated that the Compton 
continuum of the 1811 keV 7 ray leads to a correction of 
less than 1% of the raw counts in each of the ROIs. 

Similarly minor effects are expected if the counting 
rate at high energies, beyond the 0.2-2.0 MeV energy range 
discussed in the present work, differs significantly between 
the two runs. Due to the limited statistics, this high 7-ray 
energy range is not considered further here. 

The '^^Cu line at 1623 keV falls within the ROI for Ea 
= 400 keV. It contributes about 6% to the total counts 
in the ROI (fig.[S|). This problem can be mitigated by ex- 
cluding the region of this line from the ROI, and thus 
accepting a somewhat reduced statistics for the 7 ray to 
be detected. 



7 Conclusions and summary 

The effects of a weak flux of energetic neutrons - in- 
duced by energetic deuterons from elastic scattering and 
the ^H(d,n)'^IIe reaction - on an underground ultra-low 
level in-beam 7-ray spectroscopy setup has been stud- 
ied experimentally at two different a-beam energies, using 
data from a HPGe detector and a silicon particle detec- 
tor. Monte Carlo simulations and tests with a weak Am-Be 



neutron source have been performed to assess the reliabil- 
ity of the results. 

Owing to the ultra-low background setting and the 
large size of the HPGe detector, a number of (n,n'7) fea- 
tures in the background of lead-shielded HPGe detectors 
at Ej = 1.4-2.2 MeV has been described experimentally 
for the first time. 

The reliability of a GEANT4-based simulation for neu- 
tron-induced effects in a high-purity germanium detector 
has been studied, and remaining differences with the data 
have been discussed. As a result, the background in the 
planned ^H(Q;,7)^Li experiment is now well understood. 

The effects of the background given by the presently 
studied ultra-low neutron flux on the future ^H(a,7)^Li 
experiment have been evaluated, and a possible empiri- 
cal parameterization of the in-beam background has been 
performed. It was shown that this parameterization holds 
the potential to use a run at a given beam energy as a 
background monitor for a run at a different beam energy, if 
the ROIs are disjunct. Even though the background is one 
order of magnitude larger than the expected ^H(Q;,7)^Li 
signal, the present data indicate that a positive measure- 
ment may be possible. 
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